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This paper characterizes the attractor structure of synchronous and asynchronous Boolean networks induced by bi- 
threshold functions. Bi-threshold functions are generalizations of standard threshold functions and have separate 
threshold values for the transitions — > 1 (up-threshold) and 1 — >■ (down-threshold). We show that synchronous 
bi-threshold systems may, just like standard threshold systems, only have fixed points and 2-cycles as attractors. 
Asynchronous bi-threshold systems (fixed permutation update sequence), on the other hand, undergo a bifurcation. 
When the difference A of the down- and up-threshold is less than 2 they only have fixed points as limit sets. However, 
for A > 2 they may have long periodic orbits. The limiting case of A = 2 is identified using a potential function 
argument. Finally, we present a series of results on the dynamics of bi-threshold systems for families of graphs. 
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1 Introduction 



A standard Boolean threshold function tk.m ■ {0, 1}' 



{0, 1} is defined by 



1, if a{xi, . . . ,Xr, 
0, otherwise. 



> k and 



(1.1) 



where (7 (a: 1, . . . ,Xm) = \{i < i < m \ Xi = This class of functions is a common choice in modehng 
biological systems [Kauffman ('1 969| l; |Karaoz et al.| ( |20()4) l], and social behaviors (e.g., joining a strike or 
revolt, adopting a new technology or contraceptives, spread of rumors and stress, and collective action), 
see, e.g., |Granovetter| ( |TW8| ; |Bulger eraLl ( [T9g9l ); |Macy| ( [TWl) ; ^^ and Macy| ( |2007l ); |Watts| ( p002) i; 
Kempeetal.|(|2003|l]. 
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A bi-threshold function is a function ti^k^^k^^m ■ {0, 1}' 



{0, 1} defined by 



if Xi = 0, 



,m 5 



(1.2) 



Here i denotes a designated argument - later it will be the vertex or cell index. We call fc^ the up-threshold 
and the down-threshold. When fc^ = the bi-threshold function coincides with a standard threshold 
function. Note that unlike the standard threshold function in \\.\\ which is symmetric, the bi-threshold 
function is quasi-symmetric (or outer-symmetric) - with the exception of index i, it only depends on its 
arguments through their sum. 



In this paper we consider synchronous and asynchronous graph dynamical systems (GDSs), see | Mortve 
and Reidys|p007] l; paca"uley and Mortveit] ( |2509] l], of the form F : {0,1}" — > {0,1}" induced by bi- 
threshold functions. These are natural extensions of threshold GDSs and capture threshold phenomena 
exhibiting hysteresis properties. Bi-threshold systems are also prevalent in social systems where each 
individual can change back-and-forth between two states; Schelling states: "Numerous social phenom- 
ena display cyclic behavior see d Schelling] |1978 | p. 86). Among his examples is whether pick-up 
volleyball games will continue through an academic semester or die (e.g., individuals regularly choosing 
to play or not play). One can also look at public health concerns such as obesity, where an individual's 
back-and-forth decisions to diet or not — which are peer influenced, |Christakis and Fowler ( 2007| l], and 
therefore ca n be at least part i ally d escribed by thresholds — are so commonplace that it has a name: "yo- 
yo dieting" {Atkinson et al. ( 1994 1]. When fc^ > fc^, a vertex that transitions from state to state 1 is 
more likely to remain in state 1 than what would be the case in a standard threshold GDS. For the state 
transitions from 1 to the situation is analogous. This suggests that the cost to change back to state 
is great or that a change to state will occur only if the conditions that gave rise to the —J' 1 transi- 
tion significantly diminish. A company that acquires and later divests itself of a compe titor is such an 
example. Examples where k^ > k^ are commonplace. For example, | Schelling ( 1978 1] states that he 
often witnesses people who start to cross the street against traffic lights, but will return to the curb if they 
observe an insufficient number of others following behind. Overshooting, whereby a group of individuals 
take some action, and with in a short time period, a s ubset of these pull back from it, is also of interest to 
the sociology community [Bischi and Merlone ( 2009 j )] and is characterized by fc^ > fc^. 

It is convenient to introduce the quantity A ^IF- k^. The first of our main results (Theorem |3.1[) 



characterizes limit cycle structure of synchronous bi-threshold GDS (also known as as Boolean networks). 



Building on the proof for threshold functions in Goles and OUvos ( jl981 ), we prove that only fixed points 
and periodic orbits of length 2 can occur for each possible combination of fc^ and k^. Since we re-use 
parts of their proof, and also since their proof only appears in French, a condensed English translation 
is included in the appendix on page 15 The situation is very different for asynchronous bi-threshold 
GDSs where a vertex permutation is used for the update sequence. Our second main result states that 
when A < 2, only fixed points can occur as limit cycles. However, for A > 2 there are graphs for 
which arbitrary length periodic orbits can be generated. The case A = 2 is identified using a potential 
function argument and represents a (2 -parameter) bifurcation in a discrete system, a phenomenon that to 
our knowledge is novel. We also include a series of results for bi-threshold dynamics on special graph 
classes. These offer examples of asynchronous bi-threshold GDSs with long periodic orbits, and may also 
serve as building blocks in construction and modeling of bi-threshold systems with given cycle structures. 
Paper organization. We introduce necessary definitions and terminology for graph dynamical systems 
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in Section |2] The two main theorems are presented in Sections 3.1 an d|3.2[ Our collection of results on 



dynamics for graph classes like trees and cycle graphs follow in Section|4]before we conclude in Section]?] 

2 Background and Terminology 

In the following we let X denote an undirected graph with vertex set v[X] = {1, 2, . . . , n} and edge set 
e[X]. To each vertex v we assign a state E K = {0, 1} and refer to this as the vertex state. Next, we 
let n[v] denote the sequence of vertices in the 1-neighborhood of v sorted in increasing order and write 

x[v] = iXn[v](l),Xn[v](2)i ■ ■ ■ ,Xnlv]{d{v) + l}) 

for the corresponding sequence of vertex states. Here d{v) denotes the degree of v. We call x = 
{xi , X2, . . . , Xn) the system state and x[v] the restricted state. The dynamics of vertex states is governed 
by a Ust of vertex functions (/„)„ where each fy : K''-'-^^'^^ — > K maps as 

Xy{t + l)^ f,{x{t)[v]) . 

In other words, the state of vertex v at time t + 1 is given by /„ evaluated at the restricted state x[v] at 
time t. An update mechanism governs how the list of vertex functions assemble to a graph dynamical 



system map (see e.g. Mortveit and Reidys (2007 1; Macauley and Mortveit (2009 1) 



F : X" — > K" 

sending the system state at time t to that at time t + 1. 

For the update mechanism we will here use synchronous and asynchronous schemes. In the former case 
we obtain Boolean networks where 

F(a;i,...,a;„) = (/i(x[l]),...,/„(x[n])). 

This sub-class of graph dynamical systems is sometimes referred to as generalized cellular automata. In 
the latter case we will consider permutation update sequences. For this we first introduce the notion of 
X -local functions. Here the X-local function : X" — > K"^ is given by 

. . . ,a;„) = [xi,X2, ■ ■ .,fv{x[v]), . . . ,x„) . 

Using TT = (tti, . . . , 7r„) S Sx (the set of all permutations of y[X]) as an update sequence, the corre- 
sponding asynchronous (or sequential) graph dynamical system map F^r : — > is given by 

Ftt = i^TiCn) O ^7r(ri-l) O • • • O . (2.1) 

We also refer to this class of asynchronous systems as (permutation) sequential dynamical systems (SDSs). 
The X-local functions are convenient when working with the asynchronous case. In this paper we will 
consider graph dynamical systems induced by bi-threshold functions, that is, systems where each vertex 
function is given as 

■f"" ~ ■l'v.kl,ki ■~ K,kl.ki,d{v) + 1 ' 

The phase space of the GDS map F : K"" — > is the directed graph with vertex set K"^ and edge 
set { [x, F(x)) I X S if"}. A state x for which there exists a positive integer p such that F^'(a;) = x is a 
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periodic point, and the smallest such integer p is the period of x. If p = I we call x a fixed point for F. 
A state that is not periodic is a transient state. Classically, the omega-limit set of x, denoted by lo{x), is 
the accumulation points of the sequence {Y'^{x)}k>vi- In the finite case, the omega-limit set is the unique 
periodic orbit reached from x under F. 

Example 2.1 To illustrate the above concepts, take X = CirC4 as graph (shown in Figure^, and choose 
thresholds = 1 and = 3. For the synchronous case we have we have for example F(l, 0, 0, 1) = 
(0, 1, 1, 0). Using the update sequence n = (1, 2, 3, 4) we obtain Ft^{1, 0, 0, 1) = (0, 0, 1, 0). The phase 
spaces o/Ftt and F are shown in Figure^ Notice that F-^ has cycles of length 3, while the maximal cycle 
length ofF is 2. 
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Fig. 1: The graph X = Circ4 (left), and the phase spaces of (middle) and F (right) for Example |2.1| 

We remark that graph dynamical systems generalize concepts such as cellular automata and Boolean 
networks, and can describe a wide range of distributed, nonlinear phenomena. 

3 cj-Limit Set Structure of Bi-Threshold GDS 

This section contains the two main results on dynamics of synchronous and asynchronous bi-threshold 
GDSs. 



3. 1 Synchronous Bi-Threshold GDSs 

Let K — {0, 1} as before, let A — {aij) be a real-valued symmetric matrix, let (fcj)"^]^ and {k\)^^^ be 
vertex-indexed sequences of up- and down-thresholds, and define the function F = (/i ,...,/„) : X" — > 
i^" by 



1 \fxi = Q and ^ o,ijXj > kj 
if Xi = 1 and J2 ^ij^j < ^^t 



Xi otherwise. 



(3.1) 



The following theorem is a generalization of Theorem A. 1 (see appendix) to the case of bi-threshold 
functions. 
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Theorem 3.1 IfF is the synchronous GDS map over the complete graph of order n with vertex functions 
as in Equation ^3.1) , then for all x G K^, there exists s S N such that F^"'"^(a;) = F'*(a;). 

The proof builds on the arguments of the proof from|Goles and Ohvos (1981) for standard threshold 



functions (see page 15 of the appendix). Note that we can use Lemma A. 2 in its original form, but 



for Lemma A. 3 changes are needed to adapt for bi-threshold functions. The position is marked [Cross- 



reference for bi-threshold systems] in the the proof of Lemma A. 3 on page 17 Before starting the proof 
of the theorem above, we first introduce the notion of bands and give a result on their structural properties. 
This is essential in the extension of the original result. 



Let Zi e S and assume that 7^ > 3. As in the proof of |Goles and Ohvosl ( |1981 1, we set 

supp(z,) = e {0, 1, 2, . . . , T - 1} : = 1}; 

and use their partition C — {Cq, Ci, C2, . . . , Cj,}. By the assumption 7^ > 3, we are guaranteed that 
p > 1. The bi-threshold functions require a more careful structural analysis of the elements of C than in 
the case of standard threshold functions. We say that C G C is of type ah if Zj — a and zji = b where j 
and j' are the indices immediately to the left and right of C, respectively (viewed modulo T). Here we 
write niab — mabiC) for the number of elements of C of type ab. 

We claim that moi = mio. Before we prove this, observe first that the sequence (zi(0), Zi(l), . . . ,Zi{T— 
1)) can be split into contiguous sub-sequences (bands) whose states contain only isolated Os, where the 
end points have state 1, and where bands are separated by sub-sequences of lengths > 2 whose state con- 
sist entirely of Os. By the construction of C, each element C € C must be fully contained in a single band. 
Our claim above is now a direct consequence of the following lemma: 

Lemma 3.2 A band either (i) contains no element C of type 01 or 10, or (ii) contains precisely one 
element C of type 01 and precisely one element C of type 10. 



Proof: Fix a band B and let C G C be the partition containing the initial element of B. There are now 
two possibilities. In the first case, C also contains the final element of B. Then C has type 00, and any 
other partition element contained in B is necessarily of type 11. In the second case, C terminates before 
the end of B. The configuration at the end of C must then be as 

c 

m ] 

( ...00|l_l_l_..._1101 ...l|00... ) 

o 



and C is of type 01. The element C containing the index after the last element of C either goes all the 
way to the end of B, in which case it is of type 10, or it terminates before that in which case the situation 
is as in the diagram above and C" is of type 11. By repeated application of this argument, the band B is 
eventually exhausted with an element C" of type 10. All other elements of C within B not included in the 
sequence of partitions C, C and so on, must be of type 11, and the proof is complete. □ 



Corollary 3.3 moi(C) = mio(C) 
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Proof (Theorem 3.1 1: Claim: If 7, > 3 for z^ e S then L{z,, Zj) < 0. 

We can write 



J2 Li^r. h) = E (E E (^^(^ + 1) - ^^(^ - 1))) = E * 

1=1 fe=o j=i /eCfc fc=o 



V 

ik 1 



where 



'^ik = E"^.?' E + 1) ^ ^ 1)) = E"*j'^j(''= + + 1) ^ '^o.ijZj{ik - 1) . 

i=i ieCfc j=i i=i 

We need to consider ik for the four types of partition elements. As in the original proof, note that ^'^o = 
0. 

Ck is of type 00: in this case Zi{lk ~ 1) = 0, Zi{lk) — 1, Zi{lk + 2qk + 1) = and Zi{lk + 2qk + 2) = 0, 
which is only possible if 

n n 

aijZjUk - 1) > fcj and ^ ai^ZjUk + 2qk + 1) < fcj , 
i=i i=i 

which implies that '^ik < 0. 

Cfc is of type 11: this case is completely analogous to the 00 case, and again we conclude that '1'^^ < 0. 

Ck is of type 10: here z,{lk - 1) = 1, Zi{lk) = 1, Zi{lk + 2qk + 1) = and Zi{lk + 2qk + 2) = 0. This 
implies that 

E (^ij^ji^k - 1) > fct ™d E ^jCfc + + 1) < fcJ I 

3=1 3=1 

leading to < kj - kf. 

Ck is of type 01: this case is essentially the same as the 10 case, but here ^'j^. < kf — kJ. 
Using the above four cases, we now have 

n p 

Hz,, Zj) = ^ < + moo • + mil • + mio(fcJ - k\) + moi(fct - fcJ) = , 

j=0 fe=0 

where the last equality follows by Corollary |3.3| Clearly, this leads to the same contradiction as in the 
proof of Theorem A. 1 □ 



An immediate consequence of Theorem |3.1| is the following: 
Corollary 3.4 A synchronous bi-threshold CDS may only have fixed points and 2-cycles as limit sets. 
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3.2 Asynchronous Bi-Threshold GDSs 

Theorem 3.5 Let X be a graph, let n Cz Sx and let { fv)v be bi-threshold functions all satisfying A(t;) = 
^ kl < 1. The sequential dynamical system map only has fixed points as limit sets. 

As before, the graph X is finite. Note also that the per-vertex thresholds fc^ and k^ need not be uniform 
for the graph. 



Proof: The proof uses a potential function based on a construction in Barrett et al. (2006i, but see 



also Goles-Chacc et al. ( 1985 1. For a given state x e X" we assign to each vertex the potential 

P{v,x) 



Xy 1 

d{v) + 2-kl, x^ = 



Note that the quantity d{v) + 2 — fcj is the smallest number of vertex states in the local state x[v] that must 
be zero to ensure that Xy remains in state zero. Similarly, an edge e = {v, v'} is assigned the potential 



P{e = {v,v'},x)^ 



1, Xy ^ Xy' 
0^ Xy Xy' 



For book-keeping, we let rij = ni{v; x) denote the number of vertices adjacent to v in state i for i = 0, 1 
and note that ng + ni = d. The system potential P{x) at the state x is the sum of all the vertex and all the 
edge potentials. For the theorem statement it is clearly sufficient to show that each application of a vertex 
function that leads to a change in a vertex state causes the system potential to drop. 

Consider first the case where Xy is mapped from to 1 which implies that ni > k^. Since a change in 
system potential only occurs for vertex v and edges incident with v, we may disregard the other potentials 
when determining this change. Denoting the system potential before and after the update by P and P' , 
we have P = (i + 2 — fcj + ni and P' = fc^ + no which implies that 

P' -P = kl + na-d-2 + kl-ni=ki + kl~2ni-2 
<-{kl^ki)^2^A{v)^2, 

and this is strictly negative whenever A = fc"^ — fc^ < 1. Similarly, for the transition where Xy maps from 
1 to one must have rii + 1 < k^ — 1 or rii < k^ ~ 2. In this case we have 

P' -P=[d + 2-kl + ni]-[ki+ no] = 2ni + 2 - k^ ~ kl 
<2ki-A + 2-kl-ki^ A(w) - 2 

as before, concluding the proof. □ 



3.3 Bifurcations in Asynchronous GDS 

A natural question now is what happens in the case where IS. = k^ — k^ — 2 since periodic orbits are 
no longer excluded by the arguments in the proof above. The following proposition shows that there are 
graphs and choices of and fc^, such that A = 2, for which there are periodic orbits of arbitrary length. 
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Proposition 3.6 The bi-threshold GDS map over X = Circ„ with update sequence tt — (1,2,3,..., n), 
thresholds fc^ = 1 and = 3, has cycles of length n ~ 1. 

Proof: We claim that the state a; = (0, 0, . . . , 0, 1, 0) is on an (n— l)-cycle. Straightforward computations 
give that the single 1-state is shifted one position to the left upon each application of F^r until the state 
y = (0, 1, 0, . . . , 0) is reached. The image of this state is z = (1, 0, 0, ... , 0, 0, 1) which is easily seen to 
map to X. The smallest number of iterations required to return to the original state a; is n — 1, producing 
a cycle as claimed. □ 

In other words, by taking A as a parameter, we see that the bi-threshold sequential dynamical system 
undergoes a bifurcation at A = 2. 

4 Dynamics of Bi-Threshold GDSs 

4. 1 Graph Unions 

From Proposition |3.6| we see that for X = Circ„ with threshold fc^ = 1 and fc^ = 3 at each vertex, we 



obtain an (n— 1) -cycle for the update sequence vr ~ (1,2, .. . ,n). The following proposition demonstrates 
how we can combine graphs to obtain larger cycle sizes for bi-threshold SDSs with arbitrarily nonuniform 
^k^. In particular, the result applies to the case where we combine Circ„ graphs where p = n — 1 is 
prime. 

Proposition 4.1 For i — 1,2 let Xi be a graph for which the bi-threshold GDS with update sequence tt^ 
has a cycle in phase space of length Ci. Let Ui € v[Xi], and let X be the graph obtained as the disjoint 
union of Xi and X2 plus additionally the vertex w ^ v[Xi] , v[X2] with the edges {ui , w} and {u2 , w}. 
Moreover, let all thresholds of vertices in Xi and X2 be as before, and assign threshold fc^ = 3 fo w. The 
bi-threshold SDS map over X with update sequence n = (7ri|7r2|u!) [juxtaposition] has a cycle of length 
lcm(ci, C2). 

Proof: Let vertex w have k^ = 3, so that w will never transition to state 1 from state 0. Let x = 
{xi\x2\xw) be the state over X constructed from states xi and X2 on the respective Q-cycle over Xi 
and X2 with x^ — 0. The only vertices whose connectivity, and therefore induced vertex function, are 
affected by the addition of w are ui and U2- But the state transitions for ui and U2 are unaffected because 
each is predicated on cr(a;[ui]) and (j{x[u2]), respectively, and these latter two quantities are not altered 
by the state of w because that state is fixed at by construction. Hence, the phase space of X contains a 
cycle of length lcm(ci ,02) as claimed. □ 

Thus, for fc^ = 1 and k^ = 3, there exists a circle graph and permutation tt that will produce a cycle in 
phase space of length three or greater, and multiple circle graphs can be combined to produce graphs with 
large orbit cycles without modifying the thresholds of vertices in Xi and X2. 

4.2 Trees 



Propositions 3.6 and 4.1 show how periodic orbits of length > 2 arise over graphs that contain cycles. 
This section investigates bi-threshold SDS maps where X is a tree. 

To start, we first recall the notion of K-equivalence of permutations from| Macauley and Mortveit| ( |2009| 



2008 1. Two permutations tt, tt' G Sx are n-equivalent if the corresponding induced acyclic orienta- 



tions and Ott' of X are related by a sequence of source-to-sink conversions. Here, the orientation 



Bifurcations in Boolean Networks 



1 

O- 



2 

-o- 



V 

-o- 



n/2 



n n-1 n-77+l n/2+1 

Fig. 2: The tree _ff„ used in the proof of Proposition |4.2| 

is obtained from tt by orienting each edge {v, v'} G e[X] as {v, v') if v precedes v' in tt and as (f v) 
otherwise. This is an equivalence relation, and it is shown in |Macauley and Mortveit ( 2009| l that (i) for a 



tree the number of ^-equivalence classes is k{X) = 1, and (ii) that and F^/ have the same periodic 
orbit structure (up to digraph isomorphism/topological conjugation) whenever tt and tt' are K-equivalent. 
As a result, we only need to consider a single permutation update sequence to study the possible periodic 
orbit structures of permutation SDS maps over a tree X. 

The following result shows that there can be cycles of length 3 or greater for permutation SDS over a 
tree. 

Proposition 4.2 For any integer c > 3 there is a tree X on n = Ac — 6 vertices such that bi-threshold 
permutation SDS maps over X with thresholds A;^ = 1 and fc"^ = 3 have periodic orbits of length c. 

Proof: An H-tree on n = 4/3 + 2 vertices, denoted by iJ„, has vertex set {1,2,..., n} and edge set 

{77, n-r] + l}Li {{i,i + 1}, n/2 + {i,i + 1} \ 1 < i < n/2 ~ l} , 

where r] ~ /3 + 1 and jS > 1. The graph iJ„ is illustrated in Figure |4!2l 

Set (5 — c — 2 so that n = 4(3 + 2 and rj = (5 + 1. We take X — Hn as t he g raph and assign 



thresholds (fc^, k^) — (1, 3) to all vertices. By the comment preceding Proposition 4.1 we may simply 
use TT = (1, 2, 3, . . . , n) as update sequence since all permutations give cycle equivalent maps F„. 

For the initial configuration, set the state of each vertex v in the range (n/2) + l < v < n — rj + 1 
(bottom right branch) to 1 and set all other vertex states to so that 

a;(0) = (0,0,...,0, 1, 1, . ■ . , 1, 0,0,. ..0) 

start at vertex (n/2) + 1 

The number of vertices in a contiguous vertex range with state 1 will always be tj; there may be one or 
two such groups in a system state. The image of x(0) is 

x{l) = (0, 0, . . . , 0, 1,1,. ^..,1, 0, 0, ... 0) , 

start at vertex Tj 



where now the first rj — 1 vertices are in state 0, the next 77 vertices are in state 1, and the remain- 
ing vertices — all those along the bottom arm — are in state 0, as follows. Along the top arm, vertices 1 
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through r] — 1 will remain in state because all nodes and their neighbors are in state 0. Vertex ry, the 

state of the vertex incident to the crossbar on the top arm, will change to 1 because its neighbor along 
the crossbar is in state 1. For the given permutation, then, each subsequent vertex Vi in the range rj + 1 
through n/2 will change to state 1 because Xy._^ = 1 and k'^ = 1. For the bottom arm, vertex (n/2) + 1 
will change from state 1 to state because o-{x[v(n/2)+i]) = 2 < k^. For the same reason, each vertex 
Vi in the range (n/2) + 2 to n — ry + 1 will transition to state 0. Vertices from n — rj + 2 through n will 
remain in state 0. 
The next state is 



a;(2) = (0,0,...,0, M^_^ ,0,0,...0, M^_^ ) 

start at vertex start at vertex n—?7+l 



where, for the top arm, the first rj — 2 vertices are in state 0, the next rj vertices are in state 1, and the last 
vertex on the top arm is in state 0. That is, the set of 1 's along the top arm has shifted one vertex left, 
as follows. Let the set of vertices in the top arm in state 1 (in be denoted Vi through Vertex 
Vi-i will transition — )• 1 because x^. = 1. Vertex Vi will remain in state 1 because cr{x[vi]) = 3 = k^. 
Likewise u^+i through will remain in state 1. However, Vi^^i will transition to state because 

a{x[vi+ri]) = 2 < k^. We refer to this behavior as a left-shift (the analogous shift to the right is a right- 
shift). For the bottom arm, the 77 vertices (labels (n/2) + 1 through n — rj) remain in state 0. Vertex 
n — 77 + 1 transitions to state 1 because the neighbor along the crossbar is in state 1. Subsequently, vertices 
n — 77 + 2 through n transition to state 1, in turn, according to tt. 
The next state is 

a;(3) = (0,0,...,0, M^_^_^ ,0,0,...0, M^_^^ ), 

start at vertex 77— 2 start at vertex 

where the set of t] vertices in state 1 in the top arm has shifted left, and the set of r] vertices in state 1 
in the bottom arm has shifted left. The shifting process embodied in the transition from state x{2) to 
a; (3) — where there is a group of vertices in state 1 in each of the top and bottom arms — can happen a total 
of (?7 — 2) times. The state after these (r; — 2) transitions is 

a;(r?) = (l,l,...,l,0,0,...0, 1,1,. ..,1, 0,0,. ..,0). 

^ V ' ^ V ' 

start at vertex 1 start at vertex n— 277+3 

The image of x(?7) is x(0), the initial state. There are 2 + (77 — 2) + 1 state transitions, and we have a limit 
cycle of length c = 77 + 1. □ 

Of course, the proof does not guarantee that c is the minimal periodic orbit size, nor that H„ is the 
minimal order tree with a periodic orbit of this length. Additionally, there may be multiple periodic orbits 
of length c. The following proposition expands on this in the case where c > 5: there exists a tree of 
smaller order than if„ that also admits a c-cycle, namely the F-trees. 

Proposition 4.3 For any integer c > 3 there is a tree on n = 3c — 2 vertices such that bi-threshold 
permutation SDS maps over this tree with thresholds k^ = 1 and fc^ = 3 have periodic orbits of length c 
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o n 
(1)277 

Fig. 3: The tree Yn used in the proof of Proposition |4.3| 

Proof: The proof is analogous to the case of the H-tiee. We take as the graph the y-tree on n = 3/? + 1 
vertices (see Figurejs]) with /3 > 1, which has vertex set {1, 2, . . . , n} and, setting rj ~ (3 + 1, edge set 

{{i,i + 1} I 1 < i < 27? - 2} U {{i,i + 1} \ 2r] < i < {n - 1)} U {rj^n} . 

Let c > 3 with n = 3c - 2 so that X y„ (and c = /3 + 1). We assign thresholds (fc^, k^) ^ (1, 3) to all 
vertices and use update sequence tt = (1, 2, 3, . . . , n) as before. As the initial configuration, set the states 
of the /3 vertices v in the range ry < u < 2r/ — 2 (all vertices in the upper right branch except 2ry — 1) to 1, 
and set all other vertex states to to form 



x(0) = (0, 0, . . . , 0, 1,1,.^. .,1, 0, 0, ... 0) . 

Start at vertex rj 

The image of x{0) is 

= (0,0,...,0, 1, 1, . , 1, 0,0,. ..0,1), 

stait at vertex {rj — 1) 

where now the first rj — 2 vertices are in state 0, the next f3 vertices are in state 1, and the remaining 
vertices — except for vertex n — are in state 0. In the upper two branches, the initial set of /3 nodes in 
state 1 shifts left for the same reasons described in the proof of Proposition 4.2 The last vertex, n, will 
change to 1 because it is adjacent to vertex t], which has state 1. 
The image of is 

x(2) = (0,0,...,0, 1,1,. ^..,1, 0,0,...0,1,1) , 

start at vertex {ri — 2) 

where the (3 nodes in state 1 beginning at vertex rj — 2 have shifted left and vertex n ~ 1 transitions to 1 
because vertex n is in state 1. Vertex n remains in state 1 because a{x[vn]) — 3. 

The mechanics of the last state transition (the left shift of f3 vertices and nodes transitioning to state 1 
in the lower branch) repeats itself a total of /3 — 2 times, at which point the state is 

a;(/3-l) = (0,1,...,1, 0, 0, . , 0,^ 1,1, ^..,1 ), 

Start at vertex (7] + 1) start at vertex (n — /3 + 2) 
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Fig. 4: The tree X„ used in the proof of Proposition |4.4| (arrows indicate vertex labehng order). 

where the only vertex in the lower vertical branch in state is 2?], the leaf node. 

Noting that vertex rj remains in state 1 on the next transition because a{x[v,-i]) = 3, all vertices in the 
upper right branch transition to 1. Vertex 2?7 also transitions to 1, giving 

The next state can be verified to be a;(0), thus completing the cycle. The cycle length is therefore c = 13+1 
as stated. □ 

Interestingly, there is no H-tree nor y-tree that generates a maximum orbit of size 2 for thresholds 
(fc^, k^) = (1, 3). However, so-called X-trees (defined below) admit cycles of any size c > 1. 

Proposition 4.4 For any integer c > 2 there is a tree X on n ~ Ac — 3 vertices such that bi-threshold 
permutation CDS maps over X with thresholds fc^ = 1 and k'^ = 3 have periodic orbits of length c. For 
c = 1, there is a tree X on n — 5 vertices that has periodic orbits of length 1 (fixed points). 

Proof: An X-tree on n = 4/3 + 1 vertices with /3 > 1 has vertex set {1, 2, . . . , n} and edge set as 
illustrated in Figure |4] Here -q ~ /? + 1 is the unique vertex of degree 4. Note first that for any n the 
all-zero state over X„ is a fixed point. 

We treat the case c~ 2 separately; use X = X^, tt = (1, 2, 3, 4, 5), and (fc^, k^) = (1, 3). It can easily 
be verified that a::(0) = (0, 1, 1, 0, 0) is mapped to a;(l) = (1, 1, 0, 1, 1) which in turn is mapped to x{0), 
constituting a 2-cycle. 

Fix c > 3, set n = 4c — 3 and then c — f3 + 1, take as the graph X = X„ with thresholds (/c^, k^) = 
(1, 3) for all vertices, and let tt = (1, 2, 3, ... , n). 

Define the initial configuration a;(0) by assigning the (5 vertices v with ij < v < 2ri — 2 (all vertices in 
the upper right branch except 2ry — 1) to 1 and set all other vertex states to 0, that is, 

x(0) = (0, 0, . . . , 0, 1,1,.^. .,1, 0, 0, ... 0) . 

Start at vertex ?y 

The image of x(0) is 

a;(l) = (0,0,...,0, 1,1, .^..,1, 0,0,. ..0, ), 

Stan at vertex (j; — 1) start at vertex 3r; — 2 
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where now the first ?7 — 2 vertices are in state 0, the next 13 vertices are in state 1, and the remaining 
vertices in branch 2 are in state 0. In branch 3, only the vertex neighboring vertex rj transitions to state 1, 
while all vertices in branch 4 transition to state 1 because ry is in state 1. 

State x{2) is generated by a left-shift of the /3 contiguous states that are 1 in branches 1 and 2, and by a 
left-shift of the /3 -h 1 contiguous state-1 vertices in branches 3 and 4, that is, 

a;(2) = (0,0,...,0, 1,1, .^..,1, 0,0,... 0, 1, 1, .. . , 1, 0). 

start at vertex (rj - 2) start at vertex 377-8 

From a;(l) there are ^ — 2 such transitions that result in the state 

1) = (0,1,1,..., 1,0, 0,...0, 1, 1, . . . , 1, 0,0,. ..,0). 

^ V/ ' ^ </ ' 

start at vertex 2 start at vertex 3r) — /3 

The next transition results in all vertices in branches 1 and 2 in state 1 since r] remains in state 1. The 
contiguous set of ^ -|- 1 vertices in branches 3 and 4 shift left, giving 

ar(^) = (l,l,...,l, 0,0,. ..,0). 

^ >/ ' 

start at vertex 877 

The image of x{(3) is x(0), and, since /? is the smallest positive time step with this property, we have 
established the presence of a periodic orbit of length c — (3 + 1. □ 

Finally, we consider a special class of bi-threshold SDSs on trees with A;^ = 1 and A;^ = {v) = 
d{v) + 1 for each vertex v. Note that the down-threshold for each vertex depends on its degree as indicated 
by the index v in k^{v). We show that such bi-threshold SDS maps always have fixed points. In such 
systems, the state of a vertex v switches from to 1 if it has at least one neighbor in state 1, and from 1 
to if it has at least one neighbor in state 0. This is an interesting contrast to the classes of bi-threshold 
SDSs on trees discussed above which have large limit cycles. 

Let X be a tree. We choose some arbitrary vertex r £ v[X] as its root, and partition X into levels 
Xo,Xi, . . . , Xjj with respect to r such that Xq = {r}, and for any i > 0, we let Xj+i be the set of 
vertices adjacent to vertices in set Xi, but not in the set Uj^iXj. We sometimes refer to Xj as level-i 
set. Let D be the number of levels. We can also define a parent-child relationship relative to this rooted 
tree, and denote p{v) as the parent of vertex v ^ r.ln our arguments below, we use any permutation tt of 
v[X], which consists of all the vertices in Xi before those in Xi_i for each i. Our result is based on the 
following property. 

Lemma 4.5 Consider a bi-threshold SDS F^r on a tree X with an arbitrary root r and permutation n 
as defined above where k'^ = 1 and k^{v) = d{v) + I for each vertex v. Let x be any state vector and 
x' = F„{x). For each vertex v other than the root, we have x'y = Xp(^yy 

Proof: Our proof is by induction on the levels, starting from the highest, i.e., Xjj- For the base case, 

consider a leaf v E Xo- We have four cases: .x„ = Xp(^^) = 1, x^ = 0, .Xp(^,) = 1, Xy = 1, Xpi^y-^ = and 
Xy = a;p(^) = 0. It is easy to verify that in the first two cases, we have x'^, = 1 and in the latter two cases, 
we have x'^ = 0, since vertex v is updated before p{v) in tt. Therefore, the statement of the lemma holds 
in the base case for all vertices v G X^. 
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Next, consider a vertex v in some level Xj, j < D. If i; is a leaf in Xj, the lemma follows by exactly 
the same argument as in the base case. Therefore, consider the case v is not a leaf. Let wi, . . . ,Wc denote 
its children. Since level j + I vertices are updated before those in level j in tt, by induction, we have 
x'^. = Xy for each Wi. Again, we have a case similar to the base case: when vertex v is updated, it has the 
same values as its children, and therefore, takes on the state of p{v). Thus, the lemma follows. □ 

This property immediately gives us the following: 

Corollary 4.6 Let X be a tree. Let tt G Sx and let {fv)v be bi-threshold functions satisfying fc^ = 1 and 
k^iy) — d{y) + 1/or each vertex v. Any SDS map only has fixed points as limit sets. 



Proof: Without loss of generality, we take tt to be the permutation in Lemma 4.5 By applying Lemma 4.5 
it is easy to verify that for any state vector x, all the vertices in levels and 1 have the same state value in 
F{x), namely Xr- By induction on i, it is easy to verify that for any i > 1, all vertices in levels 0, . . . ,i 
have the same state value (of Xr) in F'^{x). The statement follows since all permutations for a tree give 
cycle equivalent SDS maps. □ 



5 Summary and Conclusion 

This paper has analyzed the structure of w-limit sets of bi-threshold GDS. Unlike the synchronous case, 
bi-threshold SDS maps can have long periodic orbits, and this is characterized in terms of the difference of 
the up- and down-thresholds. We also analyzed certain classes of trees. The following is a list of questions 
and conjectures for possible further research. 

5. 1 Embedding and Inheritance of Dynamics 

A fundamental question in the study of GDSs is the following: if a graph X has a graph X' as an induced 
subgraph, what are the relations between the dynamics over the two graphs? Here one has to assume that 
the vertex function, and update sequences if applicable, are appropriately related. For example, is there a 
projection from the phase space of the GDS over X to the one over X'l 

In initial computational experiments we studied the dynamics for bi-threshold GDS over trees obtained 
from, e.g. iJ-trees by adding a collection of edges - results indicate that there are several classes of 
outcomes. While this is hardly a surprise, there are clear patterns in how edges are added and the dynamics 
that result. For example, some classes of edge additions give trees that have long periodic orbits just as 
in the case of H-trees. For other classes of edge additions, however, the addition of even a single edge 
causes all periodic orbits of size > 2 to disappear Further insight into the mechanisms involved could 
shed light on the the fundamental question above. 

5.2 Minimality of Trees with Given Periodic Orbit Sizes 

Our results above on the existence of trees admitting bi-threshold SDS with given periodic orbit sizes are 
not necessarily minimal. For a given c > 1 there is an X-tree with a periodic orbit of length c, but there 
may be a smaller tree (or graph in general) which admits periodic orbits of size c as well. While we have 
obtained some insight on this via sampling, no firm results have been established. 

Note. For all computational experiments involving dynamics of SDS maps over graphs in this paper we 
used a variant of InterSim ( |Kuhlman et al.|[2(311| l. 
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A Limit Cycle Structure for Standard Threshold Cellular Automata 

This appendix section contains a condensed version of the proof from Goles and 01ivos| ( |1981| l for stan- 
dard threshold functions. We have incorporated their proof for two reasons. First, only a portion of the 
original proof needs to be adapted to cover bi-threshold systems, and in this way the paper becomes self- 
contained. Second, the original proof only appears in French, and we here provide an English version. 

Let K = {0, 1}, let A = (aij)"^^;^ be a real symmetric matrix, let 9 = (^i, . . . , On) e M", and let 
F (/i, ...,/„) : K"^ — > K" be the function defined coordinate- wise by 



Theorem A.l For all x e there exists s € N such that F'*+^(a:) ~ F'*(x). 

The proof of this theorem is based on two lemmas which are given below. Note first that since if" is 
finite, for each x E K"^ there exist s, T e N (they will generally depend on x) with T > such that 



for all < r < T. Here s is the transient length of the state x. Next define the n x T matrix X{x, T) = 



where F*(a;) = (;2i(0), . . . , z„(0)) and F"+^-i(a;) = (zi(r- 1), . . . , z„(T- 1)). In other words, z 
denotes the first periodic point reached from x (after s steps) and its period is T. The columns of X{x, T) 
are the T successive periodic points of the cycle containing z. 
In general we have 




(A.l) 



Y'+'^ix) = -F'{x) and F''+'^(a;) ^ F"(a;) 



(F^(a;),/...,F^+^-i(x))by 




¥'+\x) = (zi(Z), . . . , z„(0) for < ; < T - 1 . 



Since 



¥^{x) 



F^+^(x)=F(zi(T-l),... 



Zn{T-l)) 
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wehave^i(O) = fi{zi{T - 1), ZniT - 1)), and from F^+^+'^{x) = F(F«+'(a;)) we have 

Zi{l + 1) = Mziil), . . . , z„(0) for Z = 0, . . . ,T - 2. 

We will call Zi the i'^ row of the matrix X{x, T) and let 7^ denote the smallest divisor of T such that 
Zi{l + 7i) = Zi{l) for I a {0, . . . , T — 1}, and will say that 7; is the period of the component Zi. Clearly, 
we have Zi(l + T) = Zi{l) fori e {1, 2, . . . , n} and alH e {0, . . . ,r - 1}. Let 5 = {2:1, . . . ,2„} be the 
set of rows of X(a;, T). We define the operator L : S x 5 — > M by 

T-1 

L{zi, Zj) = aij ^{zj{l + 1) - Zj{l - l))zi{l) , 
1=0 

with indices taken modulo T. 

Lemma A.2 The operator L has the following properties: 

(i) L{zi, Zj) + L{zj, Zi) = O for i. j G {1, . . . , n} (anti-symmetry), 
(ii) Ifji < 2 then L{zi, zj) = for j G {1, . . . , n}. 

Proof: For (i), since aij = aji, we have 

T-1 

L{zi,Zj) + L{zj,Zi) = aij ^{[zi{l)zj{l + 1) - Zi{l - l)zj{l)] 

1=0 

+ [Zi{l + l)Zj{l)-Zi{l)Zj{l-l)]) , 

which clearly evaluates to zero due to periodicity. For part (ii), if 7j = 1 then the row Zi is constant and 
L{zi, Zj) = 0. If 7j = 2 then the value of Zi alternates as 

Zi{0), Zi{l), Zi{0), Zi{l), Zi{0), Zi{l) 

across the i* row, and the terms in L{zi, Zj) cancel in pairs. □ 
Let Zi G S and suppose in the following that ji > 3. We set 

suppizi) = {l€{0,...,T-l}:Zi{l) = l}, 

and write 1{1) = {1,1 + 2,1 + A, . . . ,1 - A,l - 2}. Next, set 

J 0, if there is no Zq G {0, • ■ • , T — 1} such thatX(lo) C supp(zj) 
\X{lo), otherwise. 

We define Ci as the set 

Ci = {^1 + 2s e supp{zi) ■.s = 0,l,...,qi}, 

where li is the smallest index not in Co satisfying Zi{li — 2) = and qi satisfies Zi{li + 2qi + 2) = 0. 
For fc > 2 we define the sets by 

Cfe = {^fe + 2s e supp{zi) : s = 0,l,...,qk} , 
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where Ik = h-i + r (mod T) ^ {li, . . . , Ik-i} is the smallest index for which Zi{lk — 2) = and qk 
satisfies Zi{lk + 2qk + 2) = 0. 

Since 7i > 3 (assumption), there always exists li E supp(2i) for which Zi{li ~ 2) — 0. This allows 
us to build the collection of sets C = {Cq, . . . , Cp}. By construction, C is a partition of supp(zi). The 
following lemma provides the final piece needed in the proof of the main result. 

Lemma A.3 For Zi £ S and with 7.^ > 3 we have 

n 

J2Hz„Zj) <0. 

3=1 



Proof: Using the partition C of supp(zi), we have 

n n 

n p p 71 

= E E E (^^-(^ + 1) - ^^(^ - 1)) - E E E (^^-(^ + 1) - ^^(^ - 1)) 

p 

= E ' 

where we have introduced 

n 

^-^fc = E E (^J- + - 1)) ■ (A-2) 

3=1 leCk 

If Co = then \l>io = 0, and if Co = {^o, Iq + 2, . . . Jq - 2} we. have 

^(z,G + l)-z,(/-l))=0. 

leCo 

In other words, we always have vl'jo = 0, so we assume A: > in the following. From the assumption that 
7i > 3, there exists Ck such that Ck = {hjk + 2, . . . ,lk + 2qk}, so we can re-write 'i'ik as 

= E "^J" E^'^J^^'^ + 2S + 1) - 2:j(?fc +2S- 1)) 
n n 

= ^ aijZj{lk + 2gfe + 1) - ^ aijZj{lk - 1) . 

[Cross-reference for bi-threshold systems] By the construction of Ck, we have Zi{lk + 2qk + 2) = 

and Zi{lk) = 1 which, by the definition of / in ( |A.l| i, is only possible if 



^ a,jZj{lk + 2qk + 1) < 0^, and ^ aijZj{lk - I) > 9i ■ (A.3) 
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This implies that "^tk < and we conclude that 



^L{zi,zj) = ^*ife < 



k=l 



as required. 



□ 



Proof of Theorem A.l : From Lemma A. 2 we have that L is anti-symmetric so 



^^L(z„z,) = 0. 

i=i i=i 

However, if we assume that T > 3, then there is Zj with 7i > 3 and Coroll ary 1 3 . 3 1 produce s the desired 
contradiction. We conclude that T < 2. □ 
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